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Abstract: We study the forward-backward asymmetry (Afb) an d the differential 
branching ratio (DBR) in B — > K*fi + fi~ in the presence of new physics (NP) with 
different Lorentz structures. We consider NP contributions from vector-axial vector 
(VA), scalar-pseudoscalar (SP), and tensor (T) operators, as well as their combina- 
tions. We calculate the effects of these new Lorentz structures in the low-g 2 and 
high-g 2 regions, and explain their features through analytic approximations. We find 
two mechanisms that can give a significant deviation from the standard-model pre- 
dictions, in the direction indicated by the recent measurement of Afb by the Belle 
experiment. They involve the addition of the following NP operators: (i) VA, or (ii) 
a combination of SP and T (slightly better than T alone). These two mechanisms 
can be distinguished through measurements of DBR in B — > K*{i + fi~ and Ape m 
B — > K[i + [i~. 
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1. Introduction 

To date, the standard model (SM) has been enormously successful in explaining the 
measurements of particle-physics experiments. However, recently some discrepancies 
with the predictions of the SM have been observed in B decays. Some examples 
are: (i) the values of the B®-B® mixing phase sin 2/3 obtained from different penguin- 
dominated b — > s channels tend to be systematically smaller than that obtained 
from B® — > J/ipK s [0, (ii) in B — > ttK decays, it is difficult to account for all the 
experimental measurements within the SM [[J, (iii) the measurement of the B° s -B® 
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mixing phase by the CDF and DO collaborations deviates from the SM prediction 
||, (iv) the isospin asymmetry between the neutral and the charged decay modes of 
the B — > K*l + l~ decay also differs from the SM |3J]. Though the disagreements are 
only at the level of ~ 2-3cr, and hence not statistically significant, they are intriguing 
since they all appear in b — > s transitions. 

Recently, one such discrepancy has been observed in the lepton forward-backward 
asymmetry (A FB ) in the exclusive decay B — > K*ji + n~ [^, ^J. This is especially in- 
teresting since it is a CP-conserving process, whereas most of the other discrepancies 
involve CP violation. The deviation from the SM can be seen in the differential App, as 
a function of the dilepton invariant mass q 2 . In the high-g 2 region (q 2 > 14.4 GeV 2 ), 
the A FB (q 2 ) measurements tend to be larger than the SM expectations, although 
both show similar trends. The anomaly is more striking at low q 2 (1 GeV 2 < q 2 < 6 
GeV 2 ). In the first half of this region (q 2 < 3 GeV 2 ), the SM prediction is firmly 
negative [[?J, whereas the data favor positive values. Moreover, the SM predicts a 
zero crossing in Ap^(q 2 ) whose position is well-determined and free from hadronic 
uncertainties at leading order (LO) in a s [§], [|. The measurements, on the other 
hand, prefer positive values for App,(q 2 ) in the whole g 2 -range, suggesting that there 
might not be a zero crossing. Indeed, Belle has claimed that this disagreement shows 
a clear hint of physics beyond the SM JlO| . 



It is therefore quite natural to explore new-physics (NP) explanations of Apsiq 2 ), 
and look for the effect of this NP on other observables in the same decay [11]. 
B — > K*n + n~ is described by the quark-level transition b — > sfi + fi~. This is a 
flavor-changing neutral-current (FCNC) process, and is therefore expected to play 
an important role in the search for physics beyond the SM. There have already been 
a number of theoretical studies, both within the SM [|12| and in specific NP scenarios 



13], focusing on the branching fraction and Ap-g of B — > K* jj, + /i~ . For example, 
Ref. JT4] has pointed out that Apsiq 2 ) is a sensitive probe of NP that affects the 
SM Wilson coefficients. Other observables based on the K* spin amplitudes of this 



decay are at present under active theoretical and experimental analysis [Q5l [L6l 17 



Finally, more challenging observables, such as the polarized lepton forward-backward 
asymmetry [0, have also been considered, though the measurement of this quantity 
is still lacking. 

In the coming years, the LHCb experiment [|19j will collect around 6.4k events in 
the full range of q 2 for an integrated luminosity of 2 fb _1 (a nominal one-year data 
taking). This would allow the extraction of the SM zero (if it is there) of Ap-Q with 
a precision of ±0.5 GeV 2 . Indeed a dataset of 100 pb _1 would already improve the 
world precision obtained by Babar, Belle and CDF. These measurements would also 
permit many of the additional tests for NP mentioned above. 

The decays B — > X s fi + fi~ and B — > Kji + ^ are also described by b — )■ sn + n~ , 
and hence the same new physics would be expected to affect their measurements. 
The branching ratios of these decays offer significant constraints on NP contribu- 
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tions from all Lorentz structures. The possibility of a large Afb in B — > Kfi + fi 
was considered in Ref. [pOf] , where a general analysis, allowing for all possible NP 



effects, was performed. This included vector-axial vector (VA), scalar-pseudoscalar 
(SP), and tensor (T) operators. It was shown that A FB (q 2 ) in this decay cannot 
be enhanced significantly with new VA operators, while T operators can increase 
t4fb(<? 2 ) efficiently, especially when combined with the SP new physics. 

In this paper, we apply the method of Ref. [2U| to the decay B — > K* [i + . 



That is, we perform a general analysis of NP effects without restricting ourselves 
to a specific model. Our aim here is not to obtain precise predictions, but rather 
to obtain an understanding of how the NP affects the observables, and to establish 
which Lorentz structure(s) can accommodate the observed Afb(<? 2 ) anomaly. The 
impact of NP in t4 F b(<? 2 ) may be partly washed out by integrating over q 2 , so we 
study the differential A FB (q 2 ) in the entire q 2 region. 

We find that, after taking into account the constraints from relevant measure- 
ments, there are two NP Lorentz structures that can give predictions closer to the 
low-g 2 A-fb data than the SM. The first is the case in which one adds new VA oper- 
ators. Here, the values of Afb{q 2 ) can be always positive, and hence there is no zero 
crossing. In the second, NP T operators are present, which can shift the crossing 
point to much lower q 2 values. The addition of SP operators to the T operators 
allows the results to be somewhat closer to the data. We also point out the effects 
of viable NP scenarios on the differential branching fraction dB/dq 2 . 

In section ^ we review the decay B — > K* [i + fi~ within the SM. We introduce 
new physics in section ^ by adding all possible NP operators to the effective Hamil- 
tonian. We also calculate the constraints on the coefficients of these operators, and 
present the theoretical expressions for Afb(<? 2 ) and dB/dq 2 for B — > K*ji + n~ . Sec- 
tion f| contains our numerical results for A-pB^q 2 ) and dB/dq 2 with the addition of 
specific viable NP operators. In section |5|, we summarize our findings and discuss 
their implications. Some of the more complicated algebraic expressions can be found 
in the appendix [A] 

2. B K*n + n~: Standard Model 

Within the SM, the effective Hamiltonian for the quark-level transition b — > sfi + 
is 

H s cS M = -^^{^a^^^ + ^^^KPi + m^fe)^ 

+ C 9 ^ <JrfP L b) w + Cio ^ {*fPifi) PYSM } > (2-1) 

where Pl,r = (1 =F 7s)/2. The operators Cj {i = 1,..6) correspond to the Pi in 
Ref. pi|. The SM Wilson coefficients take the following values at the scale /i = 4.8 
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GeV in next-to-next-to- leading order |7j]: 

Cf = -0.304 , Cf = 4.211 + Y(q 2 ) , C 10 = -4.103 , (2.2) 

where Cf = C 7 - C 3 /3 - 4C 4 /9 - 20C 5 /3 - 80C 6 /9, q» is the sum of the 4-momenta 
of the fi + and fi~, and the function Y(q 2 ) is given by || 

Y(q 2 ) = h(q 2 , m c ) (to x + C 2 + 6 C 3 + 60 C t 

1/4 64 
- -% 2 , m 6 ) 7C 3 + -C 4 + 76C 5 + -C 6 ) (2.3) 



Here 



4 64 \ 4 64 64 

- 0) ( C 3 + -C 4 + 16 C 5 + y C 6 J + ~C 3 + yC 5 + -C 6 . 



1 

, arctan - x > 1 

~ 2 \ 4 , . n r a/x^T 



4 / 2 \ 4 

m q ) = -- I In — - - - x 1 -- (2+ar) V |x - 1| 



. 1 + y/1 — x in 

In = x < 1 

2 - 

(2.4) 

with x = 4m 2 jq 2 . A tiny weak phase has been neglected. 

The decay amplitude for B(pi) — > K*(p 2 , e) n + {j>+) is 

(K*(p 2 , e) |^(1 - 7s)6| 5(Pi)){cf «(p_)7 M i;(p + ) + C 10 w(p_) 7 ^(P4 
-2^m b (K\p 2 ,e)\sia^(l + l5 )b\B( Pl ))u(p^v(p + )\ , (2.5) 

where we have neglected the strange-quark mass m s . The expressions for the matrix 
elements as a function of form factors are given in Ref. p2]| , and are reproduced in 
Appendix [A] for the sake of completeness. 

The double differential decay rate is given by 

d 2 T 1 2wv / A,_ l2 , v 

\M\ 2 , (2.6) 



dq 2 dcos6 2mB (8ir) 

where v = a/1 — Amf/q 2 . Here A = 1 + f 2 + z 2 — 2(r + z) — 2rz, with f = m 2 K */m 2 B 
and 2! = q 2 /m 2 B . The forward-backward asymmetry for the muons is defined by 

a cost/ - - - — / dcos9- 



dq 2 d cos 7_ x dq 2 d cos 6 



A FB (q 2 ) = -^j— ^ ^ 1 ^ (27) 

dcos9 , - , + / a cost' 



o 



dq 2 d cos 9 J _ x dq 2 d cos 9 
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Figure 1: The SM prediction for Ap^iQ 2 ) in B — > K*fi + fi and the experimental mea- 
surements from Belle. This figure is taken from Ref. |To| l. 



where 9 is the angle between the momenta of the B and the /i + in the dimuon 
center-of-mass frame. 

In Fig. ||, we show the SM prediction for 74fb(<? 2 ), along with the experimental 
measurements from Belle. From this figure, we see that the discrepancy with the 
SM is the strongest in the low-g 2 region, where the SM predicts negative values of 
^4fb(<? 2 ), as well as a zero crossing. The zero of v4fb(<? 2 ) is particularly clean, because 
at this point the form-factor dependence cancels at LO, and a relation between the 
short-distance coefficients is predicted 



He(CT(9S)) 



2m B m b ^ 



4 
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where % is the point where t4fb(<?o) = 0. Next-to-leading-order (NLO) contributions 
shift the position of this zero to a higher value: g 2 = 3.90 ± 0.12 GeV 2 0. A 
substantial deviation from this zero crossing point would thus be a robust signal for 
new physics. This can occur if the NP affects Cj S and/or Cg G , or if it changes the 
above relation itself, such as by introducing new Wilson coefficients. The present 
experimental data point towards positive values of Ap^(q 2 ) in the entire q 2 region, 
thus favoring a non-crossing solution. In the following sections, we therefore look for 
sources of NP which can give rise to this feature. 



3. B — > K*fj, + fj, : New-Physics Lorentz Structures 

3.1 New-physics operators 

We now add new physics to the effective Hamiltonian for b — > sfi + fi~, so that it 
becomes 

(3.1) 



iVA 



ISP 



eff ) 
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where W^f 1 ls gi ven by Eq. Q2.1|) , while 

Uli = -^V*V tb \R v s^P L bfL^ + R A s^P L bfi^ 5 fi 
y 2n 

+ R' v sYP R b UlnH + R'a s^Pnb fil^} , (3.2) 
n s e £ = -^V t :V tb {R s sP R bflp + R P sP R bp l5 p 

+ # s sP L 6 /2yu + # P sP L 6 /2 75 At} , (3.3) 

U T e? t = -^V*V tb {c T sa, u bpa^p + zC TE sa^bpa a0 pe^ (3.4) 



y/2- 

are the new contributions. Here, Ry, Ra, R'v, R'a, B-s, Rp, R's, R'p,Ct and Cte are 
the NP couplings. For simplicity, in our numerical analysis of the forward-backward 
asymmetry and the differential branching ratio, these couplings are taken to be real. 
However, for completeness, the expressions allow for a complex-coupling analysis. 

As was done in the SM case, one can turn the expression of the effective Hamil- 
tonian for b — > sp + p~ into a decay amplitude for B(jp\) — > K*(p 2 ) p + (p + ) /i~(p_). 
This amplitude is 

M(B->jrW) = 0£-vzv»x 



(K*(p 2 , e)\*f{l - l,)b\B( Pl )) \ (Cf + R v ) u(P-hAP + ) 

+ (C w + R A ) M(p_)7 At 7 5 w(p+)| 
r (K*(p2,e)\sy»(l + l5 )b\B( Pl )) \R' v u(p^v{ P+ ) + R' A u(p-) W v(p, 



C cS 

-2-^m b (K*(p 2 ,e)\sta, u q u (l + l5 )b\B( Pl )) u(p.)-fv(p + ) 

+ (K*(p 2 ,e)\s(l + l5 )b\B( Pl )} {R s u(p„)v( P+ ) + R P u(p„) l5 v(P + )} 
+ (K*(p 2 ,eMl- l5 )b\B( Pl )) {R' s u{p.)v{p + ) + R' P u{p.) lb v{p + )} 
+ 2C T (K*(p 2 ,e)\m^b\B( Pl )) u(p_)a^v(p + ) 

+ 2iC TE e^(K*(p 2 ,e)\sa fiU b\B(p 1 )) u( P -)a a pv(p+ 



(3.5) 



where the expressions for the matrix elements [22] are reproduced in appendix [X[ 



Note that the matrix elements are functions of 7 form factors: Ao,i, 2 (q 2 ), V(q 2 ), 
3.2 Constraints on the new-physics couplings 

The constraints on the NP couplings in b — > s fi + pT are obtained mainly from 
the related decays B — > X s p + pT and B® — > p + n~ . Due to the large hadronic 
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uncertainties, the exclusive decays B — > (K,K*)fi + fi~ provide weaker constraints 
than the inclusive decay B — > X s fi + . 

Neglecting the muon and strange-quark masses, the branching ratio of B — > 
X s ji + n~ is given by 



B(B^X sf i + fi-) = B SM + B VA [\R V \ 2 + \R A \ 2 + \R' V \ 2 + \R' A \ 2 

+ B SM - VA Re [RyCg S + R* A C 10 ] + B' SM _ VA Re(R v G 



eff\ 
7 I 



(3.6) 



+ B, 



SP 



\R S \ 2 + \R P \ 2 + \R' S \ 2 + \R' P \ + B T [\C t \ 2 + MCte\ 2 } 



where 

B.SM 



B n 



dz (1 — z 



B 



VA 



B 



SM-VA 



•Ml : J 

2 By a 



16 

z 

(1-zf 



[\Cf\ 2 + C 2 W ] + 32 (1 - z) Cf Re(C; 



f) 



B' 



SM-VA 



32 Br 



dz (l-z) jl-z 2 - 
B T = 16B VA , 

ax 

dz (1 - zf , B 



2 (l-^ 2 



SP 



AB ( 



with z = q 2 /ml. The normalization constant B is 



3a 2 B{B^X c eV)\V* h V ts \ 



dz z{l -z) 2 , (3.7) 



(3.* 



32^f(m c )K(m c ) \V^ 
where m c = m c /mb- Here f(m c ) is the phase-space factor in B(B — > X c ev) p^| : 

/(m c ) = 1 - 8m 2 + 8m° - mj! - 24m^ In m c , (3.9) 
and K(m c ) is the 1-loop QCD correction factor f24 



K\m r 



1 - 



2a s (m & ) 



3tt 



- 2 -y) (l-- c ) 2 + ^ 



(3.10) 



The branching ratio of B X s fi + fi has been measured by both Belle [2l| and 
BaBar g|. In the low-g 2 (1 GeV 2 < q 2 < 6 GeV 2 ) and high-g 2 (14.4 GeV 2 < q 2 < 25 
GeV 2 ) regions, the measurements are 

(1.49 ± 0.50±g;^) x 1(T 6 , (Belle) , 
(1.8 ± 0.7 ± 0.5) x 10~ 6 , (BaBar) , (3.11) 
(1.60 ± 0.50) x 10~ 6 , (Average) . 

(0.42 ±0.12l°$) x 10- 6 , (Belle) , 
(0.50 ± 0.25ig;g?) x lO" 6 , (BaBar) , (3.12) 
(0.44 ± 0.12) x 10~ 6 , (Average) . 



B(B ^X s fi + fi ) low 



B(B -»■ X s fx + /J,- 



/high q 2 
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The SM predictions for B(B — > X s fi + fi ) in the low- and high-g 2 regions are (1.59 ± 
0.11) x 10~ 6 and (0.24 ± 0.07) x 10~ 6 , respectively |27 . 



The branching ratio of — > fi + fi in the presence of the NP operators is 



B(B S — > ji + \i 



G%o?m\Jlr Ba 




Ami 



in 



Us 



64vr 3 



Rs — R'q 



m b + m s 



v tb v t *f 



Rp — R'p 



4m 2 , 

H_ 

ml, 



2m 



m b + m a 



H t~{C\o + Ra — R'a 



m 



B s 



, (3.13) 



The SM prediction for B(B° S //+ yr) is (3.35 ± 0.32) x 10~ 9 |§. The CDF 
experiment has reported an upper bound on this branching ratio of 4.47 x 10~ 8 at 
90% C.L. |9|. 



These two decays provide complementary information about the NP operators. 
The contribution of the SP couplings to B — > X s ji + n~ is suppressed by the small 
coefficient B SP ~ 10~ 9 , as compared to B SM ~ 10 -6 . As a result, the constraints on 
the SP coefficients from this decay are rather weak. On the other hand, the coefficient 
of the tensor couplings, Bt, is an order of magnitude larger than Bgp, while the VA 
operators interfere with those of the SM (Bsm-va)- Therefore, this decay is sensitive 
mainly to the new VA and T couplings. In contrast, the main contributions to 
B® — > fi + n~ are precisely from the SP operators: there is no contribution from the 



vector couplings Ry , the axial- vector contribution proportional to R y A ' is suppressed 
by m^/mBs, and there is no tensor piece since (0\safj, u b\Bg(p)) vanishes. 

The constraints on the new VA couplings coming from B(B — > X s fi + fi~) involve 
the interference terms between the SM and the NP. When Ry and Ra are constrained 
to be real, the allowed region in the Rv-Ra parameter space therefore looks like an 
annulus, as shown in the left panel of Fig. ^ as long as no other NP couplings 
are present. When the couplings R' v and R' A are also permitted to be nonzero real 
numbers, the allowed region takes the form of an elliptical disc, as shown in the right 
panel of Fig. [| The R' VA couplings do not interfere with the SM, so their constraints 
take the form of an elliptical disc in the R' V -R' A plane. If Ry,A are not present, the 
constraints are approximately 



(') 



< 16.8 



while if Rv a are allowed, these constraints are somewhat weakened to 



R'v\ 2 + \R'a\ z < 39.7 



./ 1 2 



(3.14) 



(3.15) 



The constraints on the tensor operators also come entirely from B(B — )■ X s fi + /i~) 
and are rather tight. We find that the allowed values of the new tensor couplings are 
restricted to 

|C T | 2 + A\C TE \ 2 < 1.3 . (3.16) 
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Figure 2: Allowed parameter space in the Rv-Ra plane when R' V a couplings are absent 
(left panel) and present (right panel). All the couplings have been taken to be real. The 
"+" corresponds to the SM. 

For the SP operators, the present upper bound on B(B® — > ji + yr) gives the 
limit 

\R S - R f s \ 2 + \R P - R' P \ 2 <0.44, (3.17) 

where we have used \V ts \ = (0.0407 ± 0.0010) Q and f Bs = (0.243 ± 0.011) GeV 
[ pip . If only Rs,p or R' SP are present, this constitutes a severe constraint on the 
NP couplings. However, if both types of operators are present, these bounds can be 
evaded due to cancellations between the Rs,p and R' s P . In that case, the constraints 
on these couplings come mainly from B(B — > X s fi + /i~), and are rather weak: 

\R S \ 2 + \R P \ 2 < 45 , R' s = R s , R' P = R P . (3.18) 

However, we shall ignore such fine-tuned situations. 

3.3 Forward-backward asymmetry and the differential branching ratio 

The double differential decay rate d 2 T/dq 2 d cos 8, calculated by substituting the ma- 
trix element from Eq. (|3.5| ) into Eq. ( [2.6|) , in turn leads to the calculation of dB/dq 2 
and Afb(<7 2 )- 

The differential branching ratio is given by 



dB G 2 a 2 1 2 / - 

where t b is the lifetime of B meson. The quantity has the form 

X SP + X V A + X T + X S p-VA + X S p-T + X V a-T 



3f 



(3.19) 



(3.20) 



where the X terms are classified according to the contributions they contain coming 
from scalars-pseudoscalars, vectors-axial vector and tensor operators. Their complete 
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expressions are given in Appendix [A]. Note that the SM contribution is contained 
inside the X terms labeled by VA. Therefore whenever the new VA operators are 
absent, the X VA , X SP -va, X VA - T terms will be referred to as X SM , X SP . SM , X SM . T , 
respectively, for clarity. 

The forward-backward asymmetry can also be written in the form 



A FB (q 2 ) = 2m B ^0- Y SP + Y VA + Y T + Y SP - V a + Y SP - T + Yva-t , (3.21) 
rB L J 

with the complete expressions for the Y terms given in Appendix [A|. As in the case of 
the X terms, whenever new VA operators are absent, we refer to the Yva, Ysp-va, Yva-t 
terms as Ysm,Y SP -sm,Y S m-t, respectively. 

Most of the qualitative features of the NP impact on the above quantities can 
be easily understood if we use simplified expressions that neglect terms proportional 
to the small quantities m\ and f at appropriate places. (Note that this may not be 
valid for extremely low values of q 2 .) With this approximation, the terms in dB/dq 2 
simplify to 

X SP ^?,{\B^ + \B 2 \ 2 )m%z\, 

X VA « 2(\C\ 2 + \G\ 2 )m 2 B \ 2 + 2{\B\ 2 + \F\ 2 )(12fz + A) 

- 4Re(FG* + BC*)m 2 B {\ - z)X , 
X p ~ | Ct 1 2 (Quadratic terms in S 3 ,5 4 ,Ti) 

+ 1 Cte 1 2 (Quadratic terms in 5 3 ,5 4 , Ti) . (3.22) 

The three interference terms, Xs P -va, Xs P -t and Xva-t vanish in this approxima- 
tion, indicating that dB/dq 2 can be thought of as the simple addition of the SP, VA, 
and T contributions. 

With the same approximations, the only surviving Y terms in Afb(? 2 ) are 

Yy A ~ -Am B fzRe(A*F + B*F^ 

Y SP -t ~ m B zRe(2BlC TE + B*C T ^j x 

((2B 3 - AT x ){z - 1) + 5 4 m|A) . (3.23) 

The chiral structure of the operators ensures that all the other terms are suppressed 
by rhi — Tni/triB ~ 0.02. 

The approximate expressions in Eqs. (|3.22|) and ( p. 23 ) imply that 

• New interactions of the type only SP or only T only always increase 0, and 
hence dB/dq 2 , but do not contribute to Y. As a result, v4fb(<? 2 ) always de- 
creases in magnitude with such new physics. 

• New VA interactions, or an SP-T combination, is required in order to enhance 
^4fb(<? 2 ) significantly, or to change its sign. 
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Note that the SP-T contribution was already considered in Ref. [32] in the context 
of the inclusive decay B — > X s l + l~ . However, it was concluded that its effect was 
basically to increase the branching ratio while leaving unchanged the integrated 
A FB . For this reason, this contribution was disregarded in subsequent papers such 
as Ref. [Q. However, as we shall show here, this type of NP can in fact shift the 
differential asymmetry Aw(q 2 ) towards the Belle data. 

In order to determine the numerical values of dB/dq 2 and Ap^(q 2 ), we need 
to calculate the form factors. The theoretical predictions for A-fb{q 2 ) are rather 
uncertain in the intermediate region (7 GeV 2 < q 2 < 12 GeV 2 ) due to nearby 
charmed resonances. The predictions are relatively more robust for lower and higher 
q 2 . We therefore concentrate on calculating A FB (q 2 ) in the low-g 2 (1 GeV 2 < q 2 < 
6 GeV 2 ) and the high-g 2 (q 2 > 14.4 GeV 2 ) regions. 

3.3.1 Form factors in the low-g 2 region 

When the initial hadron contains the heavy b quark, and the final meson has a large 
energy, the hadronic form factors can be expanded in the small ratios Aqcd/^6 and 
Aqcv/E, where Aq CD is the strong interaction scale and E is the energy of the light 
meson. Neglecting corrections of 0(a s ), the 7 a-priori independent B — > K* form 
factors [see Eqs. ( |A.1|) - ([A.5|) 1 can be expressed in terms of two universal form factors 
U(q 2 ) and£ii(g 2 ) pi &5l M 51: 



Aiiq 2 ) 

Ao(q 2 ) 
V(q 2 ) 

Tiiq 2 ) 
Uq 2 ) 



2E 



K* 



m B + m K * 
m B 

m B - m K * 
rriK* 



Uq 

2,Ek* 
m B 



m B 

2\ 



Uq 2 ) , 

Uq 2 )~W) 

Uq 2 ) , 



Uq 2 



Uq 2 )-Ui 2 )- 

Here, Ek* is the energy of the K* in the B rest frame: 



(3.24) 



E 



Is 



m B 
2 



ni 



B 



(3.25) 



The g 2 -dependence of the form factors is assumed to be 



£||(S 2 ) = £||(0) 



1 — q 2 /m 2 B 



Uq 2 ) = £l(o) 



1 — q 2 /m 2 B 



(3.26) 
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/(o) 


Ci 


C 2 


M 


0.337 


0.602 


0.258 


A 2 


0.282 


1.172 


0.567 


A 


0.471 


1.505 


0.710 


V 


0.457 


1.482 


1.015 


Ti 


0.379 


1.519 


1.030 


T 2 


0.379 


0.517 


0.426 


T 3 


0.260 


1.129 


1.128 



Table 1 : Central values of the parameters of the form factors for the B — > K* transition 
[see Eq. (pD] H 



as predicted by power counting in the heavy-quark limit. In our analysis, we take 

fll(O) = 0.16 ±0.03 , fj_(0) = 0.26 ± 0.02 . (3.27) 

The previous relations get corrections of 0(a s ) || and possible A/m b contributions. 
However, for our analysis it is sufficient to stay at LO to determine which new 
couplings can induce a clear change of behavior of A-pB(q 2 ). 

3.3.2 Form factors in the high-g 2 region 

In order to estimate Afb{<! 2 ) in the high-g 2 region, we take the form factors calculated 



in the QCD sum rule approach p2| . The z (= q 2 /m 2 B ) dependence of the 7 form 
factors is given by 

f(z) = f(0)ex P ( Cl z + c 2 z 2 ). (3.28) 

The central values of the parameters /(0), c\ and c 2 for each form factor are given 
in Table |]. In order to take into account form factor uncertainties, we have used the 
maximum and minimum allowed values of the parameters /(0), c% and c 2 as given in 
[[231. 



4. A FB (q 2 ) and dB/dq 2 in the Presence of NP 

In this section, we examine the predictions for Ayb{q 2 ) and dB/dq 2 in the presence of 
NP operators. We consider different Lorentz structures of NP, as well as their com- 
binations, and examine the implications using the constraints on the new couplings 
obtained in Sec. |3~2| . In all figures, we show Afb(? 2 ) and dB/dq 2 for representative 
values of the NP couplings. The representative values have been chosen such that 
the maximum and minimum allowed values of A-pb{<1 2 ) and dB/dq 2 , as well as cases 
with interesting variations of Ayb{<1 2 )i are displayed. The same color (type) of line 
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in all four panels of a figure corresponds to the same values of NP parameters. In 
addition, for comparison, we also show the experimental data. For this numerical 
analysis, we have taken the NP couplings to be real. 



4.1 VA new-physics operators 



From the discussion following Eq. ( |3.23| ), it is expected that NP in the form of vector- 
axial vector operators may be able to enhance Ap-B(q 2 ) or change its sign. However, 
depending on whether the NP couplings are Rv,a or R' VA , the effect on Afb(<? 2 ) will 
have different features. In this section, we shall sequentially consider the scenarios in 
which (i) only Ry,A couplings are present, (ii) only R' VA couplings are present, and 
(iii) both types of couplings are allowed. 

4.1.1 Only R v , Ra couplings present 

Fig. [5] shows the results when the only NP couplings present are Ry and Ra- The 
following remarks are in order: 

• For certain values of Ry and Ra, A FB (q 2 ) can be either always positive (a 
possible solution for the Belle observation) or always negative. That is, for 
these cases there is no zero crossing point. This is easily explained because, in 
the presence of Ry and Ra, Eq. (|2.8|) becomes at LO 



Re(Cf(q 2 )) + R v = - 2 -^Cf . (4.1) 

% 

Then Ry can unbalance the contribution from C$ s , so that there is no solution, 
and consequently no zero. The effect of Ra is simply to rescale A FB . 

In general the zero crossing can be anywhere in the whole q 2 range. The crossing 
can be negetive to positive (positive crossing) or positive to negative (negative 
crossing) . 

It is possible to have a large A FF ,(q 2 ) while being consistent with the SM pre- 
diction of the differential branching ratio dB/dq 2 . This is explained by the 
different type of contributions entering the X and Y terms in Eqs. ( p. 22 ) and 

(PD- 



The differential branching ratio dB/dq 2 can be increased in the low-g 2 region 
by up to 50%. However, in such cases, A FF ,(q 2 ) becomes highly negative at high 
q 2 , inconsistent with the current data. This suggests that, in general, dB/dq 2 
will not be affected in this scenario. 
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Figure 3: The left (right) panels of the figure show Afb(q 2 ) and dB/dq 2 in the low-q 2 
(high-g 2 ) region, in the scenario where only Ry and Ra couplings are present. The different 
(colored) curves correspond to different choices of the Ry and Ra couplings inside their 
allowed region. For comparison, the experimental data are also displayed. 



4.1.2 Only R' v , R' A couplings present 

Fig. |] shows the results when the only NP couplings present are R' v and R' A . From 
the figure, we make the following observations: 

• For certain values of R' v and R' A , the position of the zero crossing is shifted 
significantly, but it is always a positive crossing, since A FB (q 2 ) is highly negative 
in the low-g 2 region. This behavior can be understood from Eq. ( |2.8|) , which 
in the presence of Ry and R' A becomes at LO 

Re(Cf(q 2 )) - = - 2 J^pL C f . (4.2) 

kio q 

In order to counteract the contribution from Cg S , we must have \R' v R' a /Ciq\ > 
Re(C| ff ). However, this is not allowed by the present measurement of the 
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Figure 4: The left (right) panels of the figure show -Afb(<7 2 ) and dB/dq 2 in the low-g 2 
(high-g 2 ) region, in the scenario where only R' v and R' A terms are present. 



branching ratio of B — y X s /x + /x . Hence, the zero crossing is always SM-like, 
i.e. always positive, which is not favored by the Belle data. 

It is possible to have dB/dq 2 consistent with the SM, simultaneously with a 
larger Afb^ 2 ) than the SM (up to 0.6), but only near the high-g 2 end. 

dB/dq 2 at low q 2 can be enhanced by up to a factor of 2, but then A FB (q 2 ) 
would become very small. On the other hand, dB/dq 2 at low q 2 can decrease 
by up to 50%, but this would result in a large negative value of A FB (q 2 ) in this 
region. 



4.1.3 All VA couplings present 

Fig. § shows A FB {q 2 ) and dB/dq 2 when all the VA NP couplings, Rv,Ra,R'v,R'a 
are present. The following different results are obtained depending on the choice of 
the couplings: 
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Figure 5: The left (right) panels of the figure show Afb(q 2 ) an d dB/dq 2 in the low-q 2 
(high-g 2 ) region, in the scenario where both Rv,A and R' v A terms are present. 



• For certain values of the couplings, t4fb(<? 2 ) can be either always positive or 
always negative. That is, there is no zero crossing point. 

• The zero crossing can be anywhere in the whole q 2 range. It can be either 
positive or negative. 

• Particularly interesting is the case of the top curve in A-pB(q 2 ) of Fig. |5]. Here we 
see that it is possible to have large Afb{q 2 ) at low q 2 , along with the suppression 
of dB/dq 2 in this region, as indicated by the Belle data. 

• It is possible to have dB/dq 2 consistent with the SM, simultaneously with a 
larger A-pB(q 2 ) than the SM (up to 0.6) in the whole q 2 region. 

The key point here is that, in order to reproduce the current experimental data, one 
needs both Rv,a and R' V a couplings. They change A-p^q 2 ) appropriately in the low- 
and high-g 2 regions, respectively. At present, the errors on the measurements are 
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Figure 6: The left (right) panels of the figure show ^4fb(<7 2 ) an d dB/dq 2 in the low-q 2 
(high-g 2 ) region, in the scenario where both Rs,p and R' s p terms are present. 



quite large. However, if future experiments reproduce the current central values with 
greater precision, this will put important constraints on any NP model proposed to 
explain the data. 

One NP model which contains VA operators (both Ry,A and B! VA ) involves Z'- 
mediated FCNCs. A recent analysis specifically notes that the measurement of 
Ap^(q 2 ) can be explained within this model. From the above analysis, we see that 
this is one case of a more general result. 

4.2 Only SP new-physics operators 

From the discussion following Eq. ( |3.23| ), NP involving only SP operators is ex- 
pected to decrease A FB (q 2 ). Fig. |6] shows the results when all the SP NP couplings, 
Rs, Rp, R' s , R' P are present. There we see that: 

• The SP operators have unobservably small effects on Afb(<? 2 ) and dB/dq 2 . 

• There is always a SM-like zero crossing. 
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Figure 7: The left (right) panels of the figure show Ap^iQ 2 ) and dB/dq 2 in the low-g 2 
(high-g 2 ) region, when NP is present only in the form of tensor operators. 



Contribution to Afb{<1 2 ) in this scenario can in principle come from the terms 
Y SP and Y SP _ SM in Eq. ( ft.21| ). However as as can be seen from Eq. ( |A.14|) , Y S p 
vanishes identically while Ysp-sm is ^-suppressed. In addition, the couplings Rs,p 
and R' s p are strongly constrained from the upper bound on B(B S — > For 
both of these reasons, these operators have a negligible effect on t4fb(? 2 ) and dB/dq 2 . 

4.3 Only T new-physics operators 

For the case where only tensor NP operators are added, A FB (q 2 ) is expected to be 
suppressed, as the discussion following Eq. ( |3.23|) suggests. Fig. [F] shows the results 
in this scenario. The following remarks are in order: 

• A FB (q 2 ) is in general suppressed in both the low- and high-g 2 regions, as ex- 
pected. 

• The zero crossing can be anywhere in the entire q 2 range, or it may disappear 
altogether. Whenever it is present, it is always a positive crossing like the 
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SM. This shift of zero crossing shows that the mrsuppressed Y$m-t term in 
Eq. (|A.14|) is important. In the absence of this term, the zero crossing point 



would have remained the same as the SM. 
• dB/dq 2 is enhanced. The enhancement can be significant, up to a factor of 2. 

Contributions to Apsiq 2 ) in this scenario are expected from the terms Yt and 
Y$M-T m Eq. ( p.21|) . However, as can be seen from Eq. (|A.14|) , Yt vanishes identi- 



cally, while Y S m-t is m/-suppressed. On the other hand, the term X T has no such 
suppression, and it contributes to G, resulting in an enhancement of dB/dq 2 . The 
increased value of also leads to the suppression of Apsiq 2 ) in Eq. ( |3.21|) . In some 
regions of parameter space, though, the contribution of the many terms in Y$m-t is 
no longer negligible. In such cases, the zero crossing shifts and Afb(<Z 2 ) at low q 2 can 
become positive. 

4.4 Simultaneous SP and T new-physics operators 

The discussion following Eq. ( |3.23| ) suggests that if both SP and T NP couplings are 
present simultaneously, there is the possibility that Afb(<7 2 ) is enhanced or changes 
sign. In this section, we quantitatively check if such an enhancement can take the 
Ap B (q 2 ) predictions closer to the current Belle measurements. We take the couplings 
Rs, Rp, R' s , R'p, Ct, Cte to be nonvanishing, and show the results in Fig. |[ From 
the figure, we see the following: 

• There is some parameter space of couplings where Afb(<7 2 ) is positive every- 
where, i.e. there is no zero crossing. 

• The absolute value of A.-p-g{q 2 ) cannot be enhanced above the SM, except at 
very low q 2 . Even here, the enhancement is very small. 

• dB/dq 2 is enhanced. The enhancement can be significant, up to a factor of 2. 
Since the contribution to t4fb(? 2 ) here comes from two terms, the m^-suppressed 



(but not negligible) Y S m-t and the not-rhrsuppressed Y SP - T [see Eq. ( |A. 14| )] , A FB (q 



;2\ 



is now expected to be larger than in the scenario with only T operators. Though this 
trend is observed in general, the severe restrictions on the SP couplings do not allow 
A FB (q 2 ) to become significantly more than the SM in magnitude. Still, A FB (q 2 ) can 
be influenced enough to cause a vanishing of zero crossing and positive A FB (q 2 ) at 
low q 2 . 

Certain NP models can contribute to the SP-T mechanism, though there are 
caveats. In the MSSM, tensor operators in b — > s fi + fi~ are induced from photino and 
zino box diagrams. However, their couplings are subleading in tan /3 with respect to 
the Higgs penguins Tensor operators can also be induced in leptoquark models 
by tree- level scalar leptoquark exchange (and a Fierz transformation). However, the 
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Figure 8: The left (right) panels of the figure show Afb(q 2 ) and dB/dq 2 in the low-q 2 
(high-g 2 ) region, when NP present is in the form of SP and T operators. 



tensor couplings are suppressed by the ratio of the Higgs vacuum expectation value 
and the scalar leptoquark mass f4(| . 



4.5 Other combinations of VA, SP, and T operators 

The pattern of the effect of NP Lorentz structures on Ap^(q 2 ) and dB/dq 2 is now 
clear, so that the results from the remaining combinations of operators can be dis- 
cerned straightforwardly. The addition of VA operators allows for an enhanced 
Afb{c[ 2 ) and a moderate enhancement of dB/dq 2 . The addition of SP operators 
does not affect the results much due to the severe restrictions on the SP couplings. 
The addition of T operators tends to enhance dB/dq 2 strongly and decrease A FB (q 2 ) 
at the same time. In addition, we have found that specific combinations of operators, 
such as SP-T, collaborate, with the results approaching the experimental data. 

4.6 Other new-physics sources that may affect A FB (g 2 ) 

In addition to the two NP mechanisms which have been found to significantly affect 
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Afb{c[ 2 ) (new VA operators, or an SP-T operator combination), there are two ad- 
ditional mechanisms that can in principle lead to the same effect. We comment on 
them here. 

In the first mechanism, NP is assumed to affect the ordinary dipole operator 
O-j = s<T alS Pfjb F a p. In this case, the Wilson coefficient C 7 S will be modified. This 
will result in the shifting of the position of zero crossing, as can be seen from Eq. ( |2.8| ). 
Now, there has been no hint of NP in the radiative decays B — > Xgj, K^j, imposing 
strong constraints on \Cf\. Still, if the effect of the NP is to simply reverse the sign of 
Cy ff , then Eq. ( |2.8|) would not be fulfilled, and a positive A FB (q 2 ) would be produced 
for low q 2 . However, this solution can be ruled out at 3a from the decay rate of 



B — > X s £ + £~ jy]. This constraint can be evaded if the couplings Ry and Ra are 
also taken to be nonzero. Thus, if there is NP in O7 whose sole effect is to reverse the 
sign of Cy ff , and the NP couplings Ry and Ra are present, it is possible to reproduce 
the Afb(<? 2 ) data. In other words, a great many things have to happen perfectly 
for this mechanism to work. We consider this very unlikely, and so consider this 
mechanism less plausible. 

Another NP possibility, independent of those included in Eq. fl3.1|) , is the addition 



of the chirally-flipped operator 0' 7 = sa al3 P L b F a/3 . The impact of this on A FB (q 2 ) 



together with other observables, was studied in Ref. [17|. There it was found that 



A-pft(q 2 ) does not significantly deviate from the SM prediction if only this operator 
is introduced. We therefore exclude the possibility of NP giving rise to 0' 7 . 



5. Discussion and Summary 

Motivated by the recent Belle measurement of the forward-backward asymmetry 
A FB (q 2 ) in B — > K*fi + fi~, indicating a discrepancy with the SM, we calculate this 
quantity in the presence of new physics (NP) in the low- and high-g 2 regions. We 
perform a systematic model-independent analysis, allowing for new vector-axial vec- 
tor (VA), scalar-pseudoscalar (SP) and tensor (T) couplings. Using the constraints 
on the new couplings from other related decays, we determine how the NP affects 
A-pft(q 2 ) and the differential branching fraction dB/dq 2 . This allows us to compare 
the effects of different NP Lorentz structures. 

The present Belle data jllj hint at a positive A-pB^ 2 ) in the whole q 2 region, i.e. 
no zero crossing as predicted by the SM. Though we look for NP sources that can 
give rise to this feature, our analysis is more general. Indeed, the discrepancy with 
the SM prediction is only at the 2a level, and this could change with more precise 
measurements in the future. We therefore focus on identifying unique features of 
all the sources of NP, and the patterns of their effects on v4fb(<? 2 ) and dB/dq 2 . We 
observe that the effects on these two quantities are correlated, which could enable 
the discrimination between different NP sources with future data. 
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We show, through analytical approximations and numerical calculations, that 
two kinds of NP scenarios can give rise to a positive Afb(<? 2 ) in the whole q 2 range: 

• NP VA operators can enhance Afb{q 2 ) in the whole q 2 region and keep its value 
positive. Both Ry,A and R' VA couplings are necessary. The terms involving Rv,a 
can make Ap-^lq 2 ) positive at low q 2 , while the terms involving R' VA can increase 
Apsiq 2 ) above its SM value in the high-g 2 region. It is therefore possible to 
very closely reproduce the Belle data. However, in general this can also lead 
to a significant suppression of dB/dq 2 . This is because the VA operators can 
interfere with the SM operators without rhi suppression, and a destructive 
interference in dB/dq 2 would tend to enhance t4fb(<? 2 )- Still, t4fb(<? 2 ) values 
close to the Belle data and dB/dq 2 consistent with the SM predictions are also 
possible in this scenario. 

• The T operators can influence Apsiq 2 ) in the low-g 2 region sufficiently to change 
its sign and make it positive. They still cannot enhance the magnitude of the 
asymmetry significantly, since the interference of these operators with the SM 
is m^-suppressed. Moreover, the addition of these operators can only enhance 
dB/dq 2 . The simultaneous presence of SP operators allows interference terms 
between SP and T operators that are not rhrsuppressed, and tends to take 
Apsiq 2 ) closer to the Belle data. However this improvement is marginal, since 
the SP couplings are highly constrained from the upper bound on B(B S — > 

If the Belle anomaly remains in future measurements, the NP source has to be one 
of the above two (or the less plausible mechanism involving a conspiracy between Oj, 
Rv and R A operators to flip the sign of O7). We will be able to distinguish between 
them through the following observations: 

• Apsiq 2 ) at high q 2 : the scenarios with only T operators or an SP-T combination 
cannot give rise to an enhanced Apsiq 2 ) at high q 2 . Such a situation necessarily 
requires VA operators, in particular involving the couplings R' VA - 

• Correlation with dB/dq 2 : the T-only or SP-T scenarios cannot in general sup- 
press dB/dq 2 much below its SM value, while the VA scenario will be able to 
do this. 

• Correlation with A FB (q 2 ) in B —¥ Kji + ji^: within the SM, A FB (q 2 ) in B — > 
Kfi + fi~ is consistent with zero since the hadronic matrix element for the B — > 
K transition does not get any contribution from axial- vector current. For the 
same reason, new VA operators cannot contribute to A^iq 2 ) in B — >■ Kfi + ^i~ . 
Hence, the only possible contribution to the asymmetry can be from SP or T 
operators. In Ref. ||20| , it was shown that, if the NP is only in the form of SP 
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or T operators, then the additional contribution to Afb{q 2 ) is proportional to 
the lepton mass and hence is highly suppressed. However, if both SP and T 
operators are present simultaneously, then Ap B (q 2 ) in B — > Kfi + fi~ can be as 
large as 15%. Thus, if the SP-T scenario is responsible for A FB (q 2 ) in B — > 
K*fi + fi~, we will also see a large A FB (q 2 ) in B — > Kfi + fi~. This measurement 
may also help distinguishing between the T-only and SP-T scenarios. 

Even if the Belle anomaly does not persist, or shows some other features, our 
analysis enables us to recognize patterns of the effects of VA, SP and T operators on 
A FB (q 2 ) and dB/dq 2 : 

• New VA operators can interfere with the SM operators without rhi suppression, 
and hence are the only ones that can reduce dB/dq 2 substantially below the 
SM prediction. They can also interfere among themselves and with the SM 
operators to give rise to a large magnitude for Afb(<? 2 ), with either sign. They 
can influence the zero crossing point go by changing the SM relation between 
Cg S and Cf 1 that determines its value at LO. 

• SP operators always tend to enhance dB/dq 2 , since their interference terms 
with the SM or VA operators are m^-suppressed. The contribution of only SP 
operators to A-fb{q 2 ) is also m^-suppressed. Moreover, the magnitudes of the SP 
couplings are severely constrained from the upper bound on B(B S — > 
Therefore, the addition of only SP operators does not significantly affect either 
dB/dq 2 or A FB (q 2 ). 

• T operators also always tend to enhance dB/dq 2 , owing to the rhi suppression 
of their interference with the SM operators. The enhancement can be rather 
strong - up to a factor of 2 - since the constraints on the couplings are relatively 
weak. If the NP is only in the form of T operators, then the contribution 
to Afb(<? 2 ) is also m/-suppressed, though the presence of many such terms 
mean that the total contribution is not insignificant. The net effect is that the 
scenario with only T operators tends to show a large dB/dq 2 enhancement and 
^4fb(<? 2 ) suppression. 

• The simultaneous presence of SP and T operators gives rise to a qualitatively 
new feature. Though the dB/dq 2 can be enhanced, one also gets an interference 
term between the SP and T operators that is not rhi- suppressed. As a result, 
a substantial effect on Afb{q 2 ) is possible, though a large effect is not possible 
due to the restriction on the SP couplings. 
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A. Analytical Calculation of Afb(q 2 ) and dB/dq 2 



The decay amplitude for B — > K*fi + /x [Eq. ( |3.5| )1 is written in terms of matrix 
elements of the quark operators. These are [O 



(K*(p K *,e)\s^(l± l5 )b\B(p B )) = Tiq ^le*-q 



A 2 (g 2 ) 

- U V""±J -r uoj^f i -p i^yjJB ^ FK" Jut ■ q -, ; 



± ie*Jm B + m K *)Ai(q ) =f + £>x*) M e* ' 1 



where 



e^A^p^T , (A.l) 

(m B + m K * ) 



A ( 2\ m B + m K * 2 m B -m K * 2 



(tf*(pjr.,e)|s<v&|fl(pB)) = le^^ -T 1 (g 2 )e* A (p B +^) ff 

+ ^^ (ri(g 2 )-T 2 ( g 2 ))e*V 

-4f T i(g 2 )- T ^(g 2 )- , 2 ^ 2 ^ W^'iPUl (A.3) 
q \ \m B m K *) / 



(K*(p K *,e)\sia^q v (l± l5 )b\B(p B )) = 2e^ Xa e*»p x K *q° T\(g 2 ) 

± i\ e* M (m| - m 2 K ,) - (p B + p#r*) M e* • g }> T 2 (g 2 



±ie - q ^~ (m 2 B -m 2 K *) f T3iqh (A ' 4) 



(K*(p K *,e)\s(l± l5 )b\B(p B )) = T e* ■ q A (q 2 ) . (A.5) 

m b 

Here we have neglected the strange-quark mass. The matrix elements are functions 
of 7 unknown form factors: A i i2 (q 2 ), V(q 2 ), Ti 2 ^(q 2 ). 

Using the above matrix elements, the decay amplitude for B — > K*fi + yT can be 
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written as 



M 



v; s v tb 



(«(p_)yu(p+)) x 



(A.6) 



| - 2Ae^ Xa e* u p x K *q rT - iBe^ + iC e* ■ q (p B + p K » )„ + iD e* ■ q q^ j 



- 2F l e ilvXa e fV p x K ,q a - iFe^ + %G e* ■ q (p B + p K *)^ + iHe*-qq^ 
+ iB 1 (u(p-)v(p+)) e* • q + iB 2 (u(pJ)^ h v(p + )) e* ■ q 

+ 8C TE (u(p_)<j^v(p + )) | - 2T 1 e*»(p B +p K *y + B 3 e*»<f - B,e* ■ gj&.g"} 
+ 2tC T e^ Xa {-2T ie * A ( Pi? +p K «Y + fl 3 e*Y - B 4 e* ■ qp x K *q a } 



with the quantities A, B,C, D, F 1 , F,G, H, that are relevant for VA interactions, 
defined as 



F 
G 

H 



2{Cf + R V - R' v ) m B (l + k^q 2 ) + 4m b C e 7 ff (l - k 2 ) » 



C = 2(C e 9 ff + R v - R[ 
D ■ 



m B (l + k) 



v) 



A 2 (q 2 ) 4m b C e 7 



■eff 



m B (l + k) 



T 2 (g 2 ) + ^f|T 3 (g 2 )" 
(1 — k ) 



4k 



m B 



l°9 + K V ~ tiy) 7- o 7—2 x ~5 



(q 2 /m 2 B) 



(C10 + i?A + ^ \ , 

(C10 + -R^ — R'j^ 



m B (l + k) 



(A.7) 
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The quantities -81,2,3,4, relevant for SP and T interactions, are defined as 



B 3 = 2(1 - k 



2\ ^i(? 2 ) — T 2 (q 2 ) 



(q 2 /m 2 B ) 

ih = -( I ,f /• - { p^T :i (q 2 )) (A.8) 



where g = {p + + p_) and /c = rriK*/m B . 

The double differential decay rate is given by 

rf 2 r f lVv'A , 



dq 2 d cos 6 2E B (8tt) 



|Afr , (A.9) 



where v = a/1 — 4m 2 /q 2 . Here A = 1 + r 2 + z 2 — 2(r + z) — 2fz, with r = m 2 K */m 2 B 
and 2; = q 2 /m 2 B . This leads to the differential branching ratio: 

% = ^^\ V »K\ 2 rn B T B V\e , (A.10) 

where t b is the lifetime of B meson. The quantity O has the form 
1 



9 3r 



Xgp + X V A + X T + X S p-VA + X S p-T + X V a-T 

where the complete expressions for the X terms are: 



(A.11) 



Xsp = 3 \B\\ 2 m 2 B zv 2 \ + 3 I-B2I 2 m 2 B z\ 

X VA = -8 \A\ 2 m%rz (v 2 - 3) A - \B\ 2 (v 2 - 3) {Ylrz + A) - |Cf (t; 2 - 3) A 2 
+ |F| 2 ^24f^ 2 - (v 2 - 3) A^ + 16 \F 1 \ 2 m 4 B rzv 2 X 

+ \G\ 2 m B \ ^ - 6(f + l)z (v 2 - 1) + 3z 2 (V - l) - (V - 3) A 

- 3 \H\ 2 m B z 2 (v 2 - 1) A + 2Re (FG*) m\\ ^ - r (v 2 - 3) + (2* + l)v 2 - 3 

+ 6Re (FH*) m 2 B z (v 2 - l) A + 6Re (GH*) m B (f - l)z (v 2 - l) A 

- 2Re (BC*) m 2 B (v 2 - 3) A(f + z - 1) 



- 26 - 



X T = I6m% \C TE \ 2 | - ABlz (2v 2 - 3) (12 fz + A) 

- AB 3 z (2v 2 - 3) A (B 4 m 2 B (r + z - 1) + 447\) 

- V92B 3 zT x [2v 2 - 3) (2f * + r - (z - l) 2 ) - 5 2 m^ (2w 2 - 3) A 2 

- 87\A ^B 4 m 2 B z (2v 2 - 3) (3f - z + 1) + 27\ (8r (u 2 - 3) + 25^ (2v 2 - 3)) 
-192^T 2 (2v 2 - 3) (3f (z + 2) - 2(z - l) 2 ) 
+ 4m| |C T | 2 ^4A ^ 2 ^2z7\ (225 3 + 5 4 m|(3f - z + 1)) 
+ 5 3 z (B 3 + 5 4 m|(f + ^ - 1)) + 4T 2 (25^ - 8f) j + 96fT 2 
+ 48^ 2 (^B 3 + 2T 1 )(B 3 rz + 27\(z(3f - 2* + 4) - 2)) + 4fTi(S 3 + 67\) 
+ 5 2 m|^ 2 A 2 

X 5P _va = 12m B m,ARe (5 2 *F) - 12m B rhi{r - l)ARe + 12m|m^ARe 

Xsp-t = 

X yA - T = 768m B m;f7\ARe (A*C T ) + 48m B m/Re (B*C TE ) x 

2B 3 (12rz + A) + B A m 2 B \{f + z-l)+T 1 (96rz + 48f - 48z 2 + 96z + 44A - 48) 

+48m|m / ARe (C*C TE ) (^2B 3 (r + z - 1) + S 4 m|A + 47\(3f - * + 1) j (A.12) 

Note that here, we do not differentiate between the contributions from the SM 
and the new VA operators. The forward-backward asymmetry can also be written 
in the form 

A FB (q 2 ) = 2m B |^ Y SP + Y VA + Y T + Y SP - V a + Y SP - T + Y VA . T j , (A.13) 
with the complete expressions for the Y terms given as 

Y SP = , 

Y VA = -Am B rzRe(^A*F + B*F^j , 
Y T = 0, 
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Ysp-va 
Ysp-t 

YvA-T 



rhi(f + z - l)Re (B*Bi) + ra|m/ARe (B{C) , 

m B zRB{2B* x C T E + B^Ct) (2B 3 (f + z - 1) + B 4 m 2 B \ + 4Ti(3r - z + 1 
2Re (F*C T ) m, (2B 3 (f + z - 1) + £ 4 m|A + Ti(44r - 4z + 4)) 
-2Re (G*C T ) m|m { (25 3 (3rz - z 2 + z + A) + S 4 m^(r - 1)A 
+4Ti (5fz + 4f - 3z 2 + 7z + 3A - 4) ) 

+2Re (H*C T ) m 2 B m lZ (2Bz(r + z - 1) + 5 4 m^A + 4Ti(3f - z + 1)) 
-64Re (F*C TE ) m|m, ^ 3 fz + 2T X (2fz + f - (z - l) 2 + A) ' 



(A.14) 
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